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We report on measurements of quantum many-body modes in ballistic wires and their de- 
pendence on Coulomb interactions, obtained from tunneling between two parallel wires 
in a GaAs/AlGaAs heterostructure while varying electron density. We observe two spin 
modes and one charge mode of the coupled wires, and map the dispersion velocities of the 
modes down to a critical density, at which spontaneous localization is observed. Theoreti- 
cal calculations of the charge velocity agree well with the data, although they also predict 
an additional charge mode that is not observed. The measured spin velocity is found to be 
smaller than theoretically predicted. 

Coulomb interactions have a profound effect on the behavior of electrons. The low energy prop- 
erties of interacting electronic systems are described by elementary excitations, which interact 
with each other only weakly. In two and three-dimensional disordered metals they are dubbed 
quasi-particles (1), as they bear a strong resemblance to free electrons (|2|), which are fermions 
carrying both charge and spin. However, the elementary excitations in one-dimensional (ID) 
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metals, known as Luttinger-liquids (001, are utterly different. Instead, each is collective, highly 
correlated and carries either spin or charge. 

We determine the dispersions of the elementary excitations in one-dimension by measuring the 
tunneling current. It, across an extended junction between two long ballistic parallel wires in a 
GaAs/AlGaAs heterostructure created by cleaved edge overgrowth (CEO) In this geometry 
tunneling conserves both energy and momentum. Each tunneling event creates an electron-hole 
pair with total momentum hk = eBd = hqB and total energy E = \eVsD\, in which 27rh is 
Plank's constant, — e is the electron charge, B is magnetic field applied perpendicular to the 
plane of the wires, d is the distance between their centers and Vsd is the voltage-bias between 
them @. 

The rate of tunneling between the wires depends on the ease of adding an electron to one 
wire and a hole to the other, determined by the electron-hole spectral function, Ak,E- For 
weak inter-wire interactions, Ak,E is a convolution of the individual particle spectral functions, 
which encode the overlap of electrons (or holes) with the many-body modes of the coupled- 
wires. Near Vsd = 0, in the limit of temperature T ^ 0, tunneling is appreciable only if 
IqbI = I^Fu ± ^flI, allowing exchange of electrons between the Fermi-points, kpi = vrnj/2, 
where rii is electron density in sub-band i, while i = U, L stands for sub-bands in the upper, 
lower wires. At finite energies, interactions broaden the peaks of the individual particle spectral 
functions, in particular giving them a distribution of momenta. In spite of this, at E = 0, Ak,E 
is sharply peaked at k = \kpjj ± kp^] for homogeneous wires Q. Thus, as long as momentum 
is conserved in the wires and in the tunnel junction, tunneling near ^5/) = is enhanced at the 
same B >0 values as without interactions: 

-B± = |/cfu ± ^flI • (1) 
ea 

For inhomogeneous wires, the Vsd = line-shape of the spectral function encodes information 
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on the low energy momentum distribution of the many-body states ® . 

Interactions become more important as the energy associated with them increases relative to ki- 
netic energy. To increase this ratio we reduce electron density in the wires by applying negative 
voltage, Vg, to a 2/xm top gate lying on the surface of the device. Figure[l]\ shows a typical low 
energy measurement of OIt/OVg as a function of Vg and B (6). The derivative is measured in 
order to pick up only the signal from the section of the device where density is controlled by the 
gate. This is done by adding a small ac component to Vg- A zero-bias anomaly is avoided 
by setting Vsd = 100/xV. This measurement, as well as all those reported here, is performed at 
0.25K. 

Figure [l^ shows the typical behavior of B±. At high values of Vg, tunneling is appreciable 
only in a narrow range around B±. As a function of Vg, B± evolve continuously, following 
the behavior of kpi^ and kpi^, allowing us to invert Eq. [T]and extract the density in each sub- 
band, plotted in Fig. [it (|5||21). In practice, each wire contains several sub-bands for most of 
the Vc-range. Tunneling is observed only between sub-bands with the same number of nodes 
(sub-band #1^#1, #2^#2 etc.). Tunneling amongst each pair of sub-bands, one sub-band in 
each wire, gives rise to a similar set of features. 

The dispersions of the modes can be determined for every density in the regime where we 
observe the B± peaks. The dispersions are traced by the singularities of Ak,E at finite energy 
and momentum, as depicted in Figs.|2]\,B. For non-interacting electrons, depicted in Fig.|2K, 
the curves resulting from tunneling either from or to a Fermi-point, produce a total of four 
curves: two shifted copies of the dispersion in each of the two wires (|5|). 
Finite interactions split the singularities of Ak,E (Fig. E^) because of two effects. The first is 
spin-charge separation, caused by intra-wire interactions, which creates two modes for each 
non-interacting mode. The second effect is mode-mixing, caused by inter-wire interactions. 
Generally the mixed modes are carried by both wires, giving rise to four independent velocities 
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©. This results in three identical copies of each of the four dispersions. In the limit of weak 
tunneling the spin modes do not couple, and as a result each dispersion branch in Fig.|2K splits 
into a spin mode and two coupled charge modes, creating four curves near Ig^l = |A;fu — 
and two sets of three curves near Ig^l = \kpu + kpi^\ (Fig.|2^). 

The tunneling current, It{Vsd, B), is proportional to Ak,E, whose singularities are peaks for 
the parameters of our experiment (7). During a scan of Vsd and B, It{Vsd^B) changes 
abruptly when the number of modes that can be excited with the available energy and momen- 
tum changes, where Ak,E is peaked. A curve along which this happens gives the dispersion of a 
mode, E{k). In particular, the slope at Vsd = gives the dispersion velocity, v = h''^dE /dk. 
For the experimentally relevant case of weak tunneling we expect ten such intercepts, as shown 
in Fig.|2^, but only four different magnitudes of slope. 

To determine the dependence of the dispersions on density we measured dlx (Vsd, B) /dVc for 
different Vg's, ranging from Vg = OV to Vq = —3.45V. A typical result from the regime where 
each wire has a single sub-band, —3.45 <Vg< —2.9V, is shown in Fig. 131!. One can see that 
the peaks that appeared in Fig. [lj\ at B± split and move with a slope that gives an apparent 
velocity: u = dr^ {dB / dVsD)vlo=Q- Accompanying these peaks are finite-size fringes (|9|). 
A total of six slopes appears in Fig. 13^, two near B_ and four near Both B_ slopes are 
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For each scan we extracted all 



discemable slopes. The results are summarized in Figs.|2j),E, where they are plotted versus 
the density of electrons in the first sub-band of the upper wire, rim ■ In the shaded area, which 
extends up to n* ^ SOfim^^, as extracted from Fig.[lJ\, only one sub-band is occupied per wire. 
The unexpected presence of six different branches of u in Figs. l2t),E wrongly suggests that 
the coupled-wires have more than four independent modes. The error lies in assuming that the 
band-filling induced by a finite Vsd is negligible (17^ . In reality, Vsd induces charge transfer 
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between the wires, which is controlled by the mutual capacitance, endowing kp^ and kpL with 
a Vs'D-dependence. Thus, the actual excitation velocity is given by: 

V = T^^, (2) 

where ± refers to the crossing point, B±, near which is extracted. The value of 7-1- depends 
on the capacitance matrix of the wires and is calculated using a simple model @. The model 
consists of two wires of radius r, separated from each other by a distance ^ r and from a 
nearby gate by a distance Dg/2 3> d. Since we apply Vsd to the upper wire, keeping the lower 
wire grounded, the energetic cost of adding charge to the wires is given, to quadratic order in ex- 
cess electron density, 5ni, by J2i {Epi^ni + e'^c~^6nf^ + ^ e^driiC^j^drij — edn^VsD, where 
i, j run over wire indices U, L, Epi = h^kpi'^ /{2m) is the Fermi-energy, c^^ = nh/ {2e^VFi) 
(m is the band mass of electrons, vpi = Tikpi/m is the Fermi velocity in wire i) and are 
elements of the inverse capacitance matrix. In the random phase approximation Q, the first 
term in this expression is kinetic energy, the second is Coulomb interaction energy. By as- 
sumption, the inverse-capacitance of each wire to the gate {c^ and Cll ) is identical: Cq} = 
(27re)"^ log [Dc/r]. The inverse-capacitance between the wires = Clu is: c^/ = (47re)^^ 
log [1 + (Dc/dY]. Here r = lOnm and d = 30nm. Dq = 70nm is the distance of the wires to 
a parallel layer of dopants. Using Dq = 500nm, the distance to the top gate, has only minor 
influence because of the log. The dimensions are roughly MBE growth parameters and were 
not adjusted. The result of applying Eq.|2lis presented in Fig.|3l Clearly the model is successful 
when each wire has only one occupied sub-band: all three fast branches collapse on a single 
curve for rim < n* and all slow branches collapse on two curves @. 

The same model for interactions, that corrects for band-filling, allows to identify the branches 
in Fig. |3l For this we turn to the Hamiltonian of the coupled wires, which consists of a free 
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electron part and an interacting part. Taking a long length approximation and bosonizing dTTt : 



i=ls=U I'^^^is ^ J ^ s,s'=n ° 

(3) 



The sums run over all occupied sub-bands and over both spin orientations. The density of 
spin orientation s in sub-band i is riis = ni/2, is the gradient of the displacement operator 
and Pis is the conjugate momentum. 

Within this model, which neglects back-scattering, the velocities of the coupled-wire modes 
are found by diagonalizing Eq.|3] using a canonical transformation. This yields spin velocities 
equal to the Fermi- velocities. For a single mode in each wire, N = 2, the two charge velocities 
are (|77]i : 



1/2 

(4) 



Here are the charge velocities in each individual wire (72'): Vd/vpi = y 1 + Ui/{2Epi), 
where Ui = e^rii/cG is the interaction energy. Physically +/— correspond to symmetric / anti- 
symmetric excitations (illustrated in Figs.|3^,C). When the wires are identical both modes are 
carried equally by both wires, but when the densities differ, as in the experiment, the symmetric 
mode is carried primarily by the more occupied wire, the lower wire, while the antisymmetric 
mode is carried primarily by the upper wire. 

We have overlaid the result of Eq.|4lon the corrected velocities in Figs.|3K,E. The fast velocities 
follow the calculated t>c_ closely for rim < n*, attesting to the validity of the model and leading 
us to associate them with the antisymmetric charge mode. The faster Vc+, on the other hand, 
is completely absent from the data. This is to be expected near i? , where tunneling creates 
interacting electron-hole pairs which propagate together, excitations that are almost completely 
anti- symmetric. On the other hand, near none of the excitations branches should be sup- 
pressed, as they are all excited by tunneling, leaving the issue of the complete absence of the 



symmetric mode unresolved. 

Turning to the slow branches in Figs. |3t),E, we find linear dependence on the bare Fermi- 
velocities, Vsi = vpi/ fi, where fu = fh = 1-25. The linearity and the fact that /j > 1, suggest 
that these modes are the spin modes. Theoretically one expects a spin velocity equal to f^? as 
long as back-scattering is small, while finite back-scattering reduces it below vp (I73ll7?ll75l 
I76ll77ll7^ . For example, one group (I75ll73t . using Monte Carlo simulations, finds that the 
expression /pert = 1/ y^l — y2kp/ {Trhvp) ®, gives an upper bound to the ratio between the vp 
and Vg. However, a plot of /pert in Fig.|3f; shows that /pert does not account for the deviation 
of fi from unity. 

Further examination of Fig.[T]\ down to the depletion of each sub-band in the upper wire reveals 
that the continuous evolution of each set of B± peaks is replaced by a series of vertical streaks, 
dubbed localization features (LFs). These occur in a small range of Vg below Vc* = — 3.45, — 
2.55, —2.20V, for each of the three upper- wire sub-bands we observe. Each set of LFs signals an 
abrupt change in the momentum- space content of the wavef unction in the depleting sub-band. 
Above Vg* finite-size fringes for B < i?_ and B > B+ accompany the _B-|- -peaks, signifying 
that the states in both wires contain only wavenumbers higher than the Fermi wavenumber and 
implying that they are extended {9j. The location of the fringes at B < B_ and B > B^ 
indicates that the potential along the non-uniform upper wire has a hump, with a typical length 
given by the period: h/ (eABi^inged) ~ 0.75/im, consistent with the barrier the surface gate 
induces. 

Below Vg*, we find that each LF fills a broad range in B, lying roughly between the extrapola- 
tions of B±. This implies that the wavefunction of the state along the upper wire is localized. 
We are thus led to conclude that the 5-streaks signify a qualitative change in the self-consistent 
potential at Vg*, which marks a transition between an extended state and a localized state. The 
localized states appear while the more occupied sub-bands are still fully conducting, in contrast 
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to a recent study of inhomogeneous wires dTPb . 

The localization transition affects transport along the upper wire. Figure shows the two- 
terminal conductance along this wire, G(yc,B), which is quantized. The stepwise decrease 
of G{Vg, B) with density is a hallmark of ballistic transport in a wire (20). We were able to 
measure G{Vg,B) simultaneously with dlr {Vg,B) /dVc, by recording both the dc-current 
along the upper wire and the ac-component of the tunneling current ®. The positions of the 
steps in G{Vg, B), whose dependance on B is very weak, are concurrent with the localization 
transitions apparent in Fig.|4^. We thus conclude that electrons in a sub-band cease to conduct 
because of localization while their density is still finite. 

The localization transition hints that bound states come into existence over the barrier induced 
by the top gate, which we use to vary the density. The possibility of this occurring has been 
addressed in the context of the 0.7-anomaly in point contacts (l2ill22lE?ll^?b . Using a variety 
of theoretical tools it was found that, when the density is low enough, a bound state may exist 
over the barrier (E5ll26ll27lli7l) . Our measurements show clear evidence for this scenario in long 
ID channels. Finally, 0.7-anomaly-like features are observed regularly in the conductance steps 
of CEO wires similar to ours (|2?b . Further work is needed to show a direct link between the 
0.7-anomaly and the observed localization features. 
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Figure 1: (A) Gray-scale plot of OIt (Vg, B) /dVc- LFs: localization features appear for Vq < 
Vq^, marked by arrows; Fringes: Finite size fringes. (B) Zoom on the trace of B± for sub-band 
#2. At B± are replaced by LFs, drawn schematically. (C) Dependence of density in sub- 
bands #l-#3 on Vg (Ui, L^: Upper, Lower wire sub-band i). Gray box marks regime with single 
occupied sub-band in each wire (—3.451^ < Vg < —2.90V), which starts below upper wire 
density n*. 
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Figure 2: (A) & (B) Illustration of the position of the singularities of Ak,E, with eVsn replacing 
E and B replacing k. (A) For non-interacting electrons there are four curves: / - copies of 
lower wire dispersion; // - copies of upper wire dispersion. (B) With interactions there are ten 
curves: three duplicates of each mixed charge mode (marked 3,4), two copies of each spin mode 
(marked 1,2). (C) Numerical derivative, with respect to B, of the measured OIt (Vsd, B) /dVa 
at Vg = —3.00V. Finite size fringes appear for B > B^ and B < B_. Marked are all extracted 
slopes, which are offset for clarity. Triangles- slopes extracted near B_, giving: 
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(empty). (D) & (E) Apparent 



velocities, u, versus density, rim. Each wire has a single occupied sub-band in the shaded region 
(riui < n*). (D) Dependence of m's calculated for small slopes on density. Overlayed are curves 
calculated by setting v to Vc- in Eq.|2land solving for u"^. (E) Dependence of m's calculated for 
large slopes on density. Overlayed are curves calculated by setting f su,l = ^'fu.l/I-SS in Eq.|2l 
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Figure 3: Excitation velocity versus density. (A) Velocities obtained from Fig. |2t)- Curves 
are the charge velocities Vc- (green) and Vc+ (black) (see Eq. El). (B)/(C) Illustration of the 
symmetric / anti-symmetric coupled- wire mode (+: excess positive charge, — : excess negative 
charge). (D) Velocities obtained from Fig. |2^. The lines are = vfu/1-25 (magenta) and 
"^^sL = vpi^/1.25 (red). The scale is the same as in (A). (E) Plot ofvp/v for the velocities in (A), 
(B) where w ^ is calculated from rim (for ffasf ^fast» ^siow ^siow & ^fast) or from the density in 
the first sub-band of the lower wire (for f and Vsl)- The dashed blue line is /pert (see text). 
The red and magenta curves from (D) overlap here. 
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Figure 4: (A) Two-terminal conductance of the upper wire as a function of Vq, measured by 
applying a dc-voltage bias of 100/iV along the wire. The step height deviates from the universal 
2e'^/h because of the indirect upper wire contact (t28\l29l . This measurement depends very 
weakly on B, which is 1.3T here. (B) Simultaneous measurement dlx {Vg, B) /dVc, showing 
that localization is concurrent with the conductance drops. 
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Supporting text 



The actual values of the charge velocity, f „ and the spin velocity, Vg, depend on microscopic 
details and are very difficult to determine, both theoretically and in experiment. Of particular 
interest is their dependence on n, which directly controls the ratio between the Coulomb inter- 
action and kinetic energy. A one-dimensional charge mode, which resembles a charge density 
wave, travels with a velocity that is strongly affected by the Coulomb interaction: Vc = vp/g. 
Here the Fermi- velocity is a measure of electron density and g is a measure the relative strength 
of the Coulomb interaction. For repulsive interactions < g < 1, while for non-interacting 
electrons g = I. Within the random phase approximation, which gives reliable estimates for g 
if backscattering is weak, it is found that as n is reduced, g decreases. 

The propagation velocity of the spin modes, Vs, is related to exchange interaction. According to 
theory, Vs is suppressed for very strong repulsive interactions, where it is difficult for adjacent 
electrons to exchange places, leading to Vg <^ vp (EZH^EHEIl- The main text cites an 
example for the suppression factor of Vg relative to t>i? in one model (1521 15*51) . which is given by 



l//pert = Y 1 ~ ^2kF/ {j^TT'Vf)- Hcrc Vfe is the Fourier transform of the two-body interaction 
potential in a single wire in the presence of a gate. It is given by: 



where Kq is a Bessel function, is the Bohr radius in GaAs and Dq is the distance to the top 
gate, 500nm. 

A remark is due on the number of expected singularity branches crossing the Vsd = axis in a 
measurement of dip {Vsd, B) /dVc- In principle, even in the absence of inter- wire interactions, 
we expect for each charge mode an extra feature with opposite slope and very small amplitude 
(unless interactions are very strong). This gives an extra copy of each charge mode branch 
near B_ . The reason is that forward scattering has a contribution from interactions between left 
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Figure SI: Schematic of the device, cleavage-plane facing front, perpendicular to B. Depicted: 
2/im-wide top gates (Gi-s), ohmic contacts to 2DEG (Oi-s), 20nm-wide upper wire at edge 
of 2DEG, 30nm-wide lower wire and 6nm insulating AlGaAs barrier, Vg'- voltage applied to 
control density, lu' current through upper wire. It' tunneling current. Also illustrated: U{x), 
gate-induced potential, and electrochemical potential in the upper wire, yU. 



and right movers, causing an electron tunneling into one branch of movers to induce a density 
modulation in the other branch as well. 

Our device consists of two parallel wires separated by a tunnel junction at the edges of two 
quantum wells in a GaAs/AlGaAs heterostructure created by cleaved edge overgrowth SS6[ . The 
band structure is such that besides the wires along the edge, the upper well is also occupied by 
a two-dimensional electron gas (2DEG), which we use to contact the wires. The measurements 
were conducted in a '^He-fridge at a base temperature of 0.25K. After the sample cools down, 
it is illuminated by infra-red light, which ionizes impurities and increases the overall electron 
density in the device. 

Tungsten top-gates, 2yum-wide, lying 500nm above the tunnel junction and 2/im apart, control 
the density in sections of the device. Both the tunnel junction and the upper wire are delimited 
by applying voltage to two peripheral gates, Gi & G3 in Fig.|STl lying lO/im apart. Bias voltage 
Vg is applied to a central gate, G2 in Fig. ISU to vary the density in the central section of the 
device. To increase sensitivity to processes affected by G2, we measure the derivative of It with 
respect to Vg- To this end we add a small ac-component to Vg (a few mV at a few Hz), apply 
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finite VsD to contact Oi and pick up the resulting current at contact O3 with a lock-in amplifier. 
An additional contact, O2, can be grounded in order to measure the current along the upper 
wire, J(7, but is left floating otherwise. To measure G{Vc, B) we apply voltage V to contact Oi 
(cf. Fig. IS II) and ground contact O2. This gives net bias Vsd = cuV relative to the lower wire 
(where < a < 1) and a bias drop of V along the upper wire. 

In the Vc-range being studied we do not observe transitions between a sub-band in one wire and 
more than one sub-band in the other wire. This hints at a selection rule, which we conjecture 
arises from the similarity of the wavefunctions in the two wires in the planes perpendicular 
to them (liSTI) . According to the rule, the overlap of wavefunctions in different wires with a 
different number of nodes perpendicular to the wires is small and suppresses the transitions 
between them. In identical wires this selection rule would be absolute, because wavefunctions 
from different sub-bands would be orthogonal. 

The finite width of the wires in the direction perpendicular to B distorts the dispersions slightly. 
Equation 1 in the main text is precise only as long as a single sub-band is occupied in each 
wire. When higher sub-bands are occupied, finite B flattens the dispersions (^1, changing the 
occupations and generally making it energetically favorable to occupy lower sub-bands at the 
expense of higher ones. We ignore this distortion as it gives an error that is comparable to the 
accuracy of the measurement. 

In the main text we describe a model that allows to calculate the actual excitation velocities 
from the measured slopes. According to the model, the value of 7± appearing in Eq. 2 in the 
main text is given by: 

7± = ^ [c±Cu (cl + Czp/2)] [(c+ + c_) (cl + Cu) + c+c_ + 4clCu]""^ , (SI) 

where c^^ = (^Cq^ ± cj}^ /2, while the rest of the quantities here are given in the main text. 
The model has a limited range of validity. Above n*, as more sub-bands are occupied, it ex- 
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aggerates the voltage-induced band-filling. In this regime, the sparsely occupied sub-bands are 
filled instead of the more populated ones, because they are more compressible. Thus in this 
regime Eq. 2 corrects the w's too much, bringing the corrected velocities too low, as can be seen 
for in Fig. 3. To obtain the fits shown in this figure we used the same values of r and 

d as for the band-filling correction. For Do we used here 500nm, the distance to the metallic 
top gate, rather than the distance to the dopant layer, because typical timescales for the reaction 
of that layer are too slow to influence the dynamics of the modes. 
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